NASA-TM-111404 


AIAA 96-1752 

/y 7/ 

B 77 7^ 


COMPUTATION OF ACOUSTIC WAVES THROUGH 
SLIDING-ZONE INTERFACES USING AN 
EULER/NAVIER-STOKES CODE 

by 

Christopher L. Rumsey 
NASA Langley Research Center 
Mail Stop 128, Hampton, VA 23681-0001 


Paper presented at the 

2nd AIAA/CEAS Aeroacoustics Conference 
May 6-8, 1996 
State College, PA 




AIAA 96-1752 
2nd A1AA/CEAS Aeroacoustics Conference 
May 6-8, 1996 / State College, PA 


COMPUTATION OF ACOUSTIC WAVES THROUGH SLIDING-ZONE 
INTERFACES USING AN EULER/NAVIER-STOKES CODE 

Christopher L. Rumsey* 

NASA Langley Research Center , Mail Stop I28 t Hampton , VA 23681-0001 


Abstract 


The effect of a patched sliding-zone interface on the transmission of acoustic waves is examined for 
two- and three-dimensional model problems. A simple but general interpolation scheme at the patched 
boundary passes acoustic waves without distortion, provided that a sufficiently small time step is 
taken. A guideline is provided for the maximum permissible time step or zone speed that gives an 
acceptable error introduced by the sliding-zone interface. 


Introduction 

The prediction and control of ducted fan noise are im- 
portant elements of the NASA Advanced Subsonic 
Technology Noise Reduction Program. Current meth- 
ods of prediction rely extensively on field measure- 
ments and analytical scaling techniques. As computers 
continue to become more powerful, Euler and Navier- 
Stokes computer codes for ducted-fan noise prediction 
have become increasingly affordable. Recent advances 
in algorithms that can enhance the efficiency of 
time-accurate Euler and Navier-Stokes computations 
also make difficult computations, such as three-dimen- 
sional (3D) rotor-stator interactions, better suited for in- 
clusion in future engine design cycles. 

Previous work by Rai 4 and Gundy-Burlet et al. 5 dem- 
onstrated the feasibility of using Navier-Stokes compu- 
tations for time-accurate rotor-stator interactions. 
Gundy-Burlet et al. utilized a combination of overlap- 
ping and patched grids; the motion of the rotor relative 
to the stator was accomplished by “sliding” the rotor 
grid system past the stator grid system and utilizing a 
non-conservative linear interpolation to transfer infor- 
mation between the two grid systems. Hall and 
Delaney 6 employed a similar patched sliding-zone inter- 
face strategy to compute ducted prop-fan flows. Chen 
and Chakravarthy 7 utilized patched sliding-zone inter- 
faces to perform rotor-stator computations. They used a 
simple piecewise-constant projection of flow variables 
between grid zones, with an area- weighting strategy. Ja- 
nus and Whitfield 8 utilized localized grid distortion to 
pass information between zones that move relative to 
one another in a prop-fan simulation. Rather than em- 


ploy interpolation, grid points near the zone interface 
were distorted, then “clicked” to new positions when ap- 
propriate. In refs. 4-8, the focus of the computations 
was the prediction of global aerodynamic characteris- 
tics. Rangwalla and Rar compared the numerically cal- 
culated tonal acoustics with theoretical values for a two- 
dimensional (2D) rotor-stator interaction. Emphasis was 
placed on the effects of boundary conditions and bound- 
ary extent; however, the effect of the patched sliding- 
zone interface on the accuracy of simulating the passage 
of acoustic waves was not explored. 

Because accurate prediction of acoustic waves is es- 
sential to any noise-prediction analysis, the effect of the 
sliding-zone interface on the passage of such waves 
must be addressed. In the present paper, the effect of a 
patched sliding-zone interface similar to that employed 
in refs. 5-7 is examined for several 2D problems, as well 
as for the passage of typical rotor-stator interaction 
modes through a 3D duct with a rotating zone. The ef- 
fects of time step and speed of the moving zone are ex- 
amined; an engineering rule of thumb for maximum 
permissible time step or zone speed is developed. 

Governing Equations 

The compressible thin-layer Navier-Stokes equations, 
written in an inertial reference frame in generalized co- 
ordinates, are 


ae . a(f-/,.) . afc-cj a( 

dt an a; 


= 0 


( 1 ) 


where 
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and Q is the vector of conserved variables given as 
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(II) 


p 


Q = 


pu 
pv 
P w 


The inviscid flux terms are 
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and 


= 




(14) 
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and t and q are the shear stress and heat flux terms, re* 
spectively. When the thin-layer assumption is invoked, 
the viscous cross-derivative terms in \ , ri , and £ are 
(6) neglected. 

For the applications here, only the inviscid (Euler) 
equations are employed. The equations are closed by the 
equation of state for a perfect gas 


The variable / represents the Jacobian of the transforma- 
tion 


P = 


(Y-l) 


[ ,22 2 . 
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e "P 2 ^ 


(15) 
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and the contravariant velocities are given by 


and 


V = q r w + r| v v + r\.w + r\ r 
W= ^n + ^v + ^ + C, 


The viscous terms are 


(7) Description of the Code 

The computer code CFL3D 10 solves the 3D time-de- 
pendent thin-layer Navier-Stokes equations with an up- 
wind finite-volume formulation. This code can solve 
flows over multiple-zone grids that are connected in a 

(8) one-to-one, patched, or overset manner, and can employ 
grid sequencing, multigrid, and local time stepping in 
accelerating convergence to steady state. Upwind-bi- 

W ased spatial differencing is used for the inviscid terms, 
and flux limiting is used to obtain smooth solutions in 
the vicinity of shock waves. Viscous terms, when em- 
ployed, are centrally differenced. The equations are 
( 10 ) solved implicitly with the use of a three-factor approxi- 
mate factorization (AF). Either the flux-difference-split- 
ting (FDS) method of Roe 11 or the flux-vector-splitting 
(FVS) method of Van Leer 12 can be employed to obtain 
fluxes at the cell faces. The FDS approach is used for all 
results in this paper. 
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Time-Advancement and Subiteration Algorithm 

The CFL3D code is advanced in time with an implicit 
AF method. The implicit spatial derivatives are first-or- 
der accurate, which results in block-tridiagonal inver- 
sions for each sweep. However, for solutions that utilize 
FDS the block-tridiagonal inversions are further simpli- 
fied with a diagonal algorithm; when the viscous terms 
are used, a spectral radius scaling is employed. 

Second-order temporal accuracy for a single-step AF 
scheme is forfeited for unsteady computations with 
these simplifications to the left-hand side. One method 
for recovering the desired accuracy is through the use of 
subiterations. Two subiteration strategies have been im- 
plemented in CFL3D. These strategies were explored in 
detail by Rumsey et al. 3 The method employed for the 
computations in this paper is termed the “pseudo time 
subiteration” (or “x -TS”) method, which uses a second- 
order-accurate temporal discretization. In the literature 
(e.g., Venkateswaran and Merkle 13 ), the t-TS method 
is also referred to as the “dual time” method because a 
pseudo time is used to iterate to the next desired physical 
time. The chief advantage to the x -TS subiteration strat- 
egy is that it frees the user from numerical constraints on 
the time step; an appropriate physical time step can be 
selected to resolve the physics of the flow in question. 

Dynamic Patched-Grid Algorithm 

A patched-grid interface implies that two zones share 
a common interface at which the grid points do not nec- 
essarily connect in a one-to-one manner. For grids in rel- 
ative motion to one another (for example, when a grid 
zone around a rotor slides past a grid zone around a sta- 
tor), dynamic patched interfacing is a relatively easy 
way to handle the transfer of time-accurate data between 
zones. 

In CFL3D, the patching is accomplished non-conser- 
vatively in the following way. First, at the patched inter- 
face, the problem is reduced by 1 degree of freedom 
because data are used only in the two planes of cell cen- 
ters in computational space nearest each interface to up- 
date the ghost-cell boundary conditions in the other 
zone. For example (as shown in figure 1 in two dimen- 
sions), if the k = kdim boundary of zone 1 interfaces with 
the k = 1 boundary of zone 2, then zone 2 obtains infor- 
mation for its two ghost cells from the k - kdim - 1 and 
k = kdim - 2 cell centers of zone 1 ; zone 1 obtains infor- 
mation for its two ghost cells from the k = 1 and k = 2 
cell centers of zone 2. 

Interpolation is accomplished dynamically at each 
time step (every time the moving zone changes posi- 
tion). For each face center point on the interface in the 
“to” zone, a corresponding real-valued index location is 
determined in the “from" zone. The primitive variables 
are then interpolated to this real-valued index location. 
A simple 2D example is given in figure 2; zone 2 acts as 
the “to” zone, and only one row of ghost cells is shown 
in the “from” zone for clarity. In this example, the center 
point of the cell j = 1 in zone 2 lines up with r| = 0.6667 
in zone 1, which lies in the cell j = 1 of zone 1. There- 
fore, this particular ghost-cell boundary condition is in- 
terpolated by using the primitive variables at = 1 in 
zone 1 and an appropriately weighted fraction of the r| - 
direction gradient at y = 1 . Further details of the patched- 
grid algorithm are given in Biedron and Thomas. 14 

When periodic boundary conditions are applied and 


Patched interface (gap shown 
for clarity only) 



center 


Fig. 1. Sample 2D patched interface. 

Zone 1 (“from”) Zone 2 (“to”) 



Fig. 2. Sample 2D patched interface, with zone 2 as “to” zone. 

one grid zone slides relative to another, copies of both 
the moving and nonmoving zones must be supplied, 
which are translated or rotated appropriately to ensure 
that each “to” zone has a “from” neighbor at all times 
and vice versa. To avoid the necessity of making too 
many copies over long times, the current algorithm also 
periodically rotates or translates the moving block and 
its solution whenever its movement exceeds a specified 
limit. 

Results 

Several 2D and 3D model problems are examined to 
study the effect of a patched sliding-zone interface on 
the passage of acoustic waves. Unless otherwise noted, 
three x -TS subiterations with three levels of multigrid 
are used for the results in this paper. 
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2D Axiallv-Moving Acoustic Waves 

The first test case involves acoustic waves that are 
moving axially in two dimensions (jc , z ). At the left end 
of a computational domain, waves are generated by per- 
turbing the free-stream pressure levels in accordance 
with 

P = p 0 + /icos [ Y _r ) cos ( (16) 

where p 0 is the nondimensional free-stream static pres- 
sure (= 1/y) and h is taken as 0.04518. Time is nondi- 
mensionalized by unit length and speed of sound. The 
free-stream is quiescent. 

This test case is patterned after a computation by 
Khan, 15 who used a perturbational form of the Euler 
equations to compute the nonlinear behavior of a plane 
acoustic wave that was propagating axially; the results 
were compared with the theoretical predictions of 
Pierce 16 at locations prior to where nonlinear steepening 
causes shock formation. In the current study, a sinusoi- 
dal transverse variation of pressure (in the z direction) is 
also imposed to test the effect of a sliding-zone interface 
(no effect is realized otherwise). Comparison with the 
theoretical results is made at a transverse location at 
which the acoustic wave front is aligned with the z-axis. 
The variation in the z direction is gradual enough to al- 
low the wave to approximate a plane wave at this loca- 
tion. 

The computation is performed by using either one or 
two zones. In the latter case, zone 2 slides in the z direc- 
tion along one end of zone 1 ; a patched sliding- zone in- 
terface connects the two zones. (See figure 3.) The upper 
and lower boundaries are treated as periodic interfaces, 
and the right boundary employs a characteristic bound- 
ary condition. To minimize the possibility of spurious 
reflections from the right boundary, the grid is stretched 
in the x direction past x = 1 .6345 . 


Sliding interface 

/ . t 

z= 3.2 


/Zone 1 

V Zone 2 

* 

a 

. _ . A 



2=0 


x=0.817 » x=3.4 

ressure perturbation at x=0 


Fig. 3. Sketch of 2D case for axially-moving acoustic waves. 

The parameter T t is taken as 0.55164 (which corre- 
sponds to a frequency of 0.6 kHz in air), and T. is taken 


4 


as 3.2. This latter parameter yields one complete period 
of pressure fluctuation between z = 0 and z = 3.2 . 

Prior to evaluating the effect of the sliding-zone inter- 
face, parametric studies were used to determine the grid 
size and time step necessary to capture the physics of 
this flow field with the CFL3D code. Figure 4 shows the 
effect of varying the number of cells in the z direction 
from 8 to 64; 16 to 32 grid points are sufficient to accu- 
rately depict the sinusoidal variation in the acoustic 
pressure. 



Z 


Fig. 4. Effect of number of points in z direction on pressures at 
x = 0.83 , with Ax = 0.01277 , At = 0.005 , and T = 3.6. 

The effect of Ax is shown in figures 5 and 6. The 
wavelength in the x direction for this case is 0.552. Us- 
ing a constant Az spacing of 0.1 (32 points per wave in 
the z direction), the Ax spacing is varied from 0.02554 
(approximately 22 points per wave) to 0.006385 (ap- 
proximately 86 points per wave). Figure 5 shows the 
pressure as a function of x/x along z = 0 , where x is 
the value of the distance at which the shock forms: 



(' indicates a dimensional quantity). In this equation, 
k' = a>'/a ' , to' is the frequency in rad/s, p' is the den- 
sity, and a' is the speed of sound. For this case, 
x = 1.6194 . Figure 6 shows the amplitude of the funda- 
mental and the first two harmonics as a function of dis- 
tance. The computational curves in this figure are 
generated by establishing a periodic solution, comput- 
ing over a nondimensional time of 1 .1, and using a Fou- 
rier decomposition on the time histories of pressure at 
each spatial location. Three higher harmonics are com- 
puted; only the first two are shown in figure 6 for clarity. 
The predicted fundamental and the first two harmonics 
agree well with the theoretical predictions of Pierce 16 
with 43 or more points per wave. With only 22 points 
per wave, significant error is realized as x/x approaches 
1 , particularly for the higher harmonics. 
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Fig. 5. Effect of number of points in jc direction on pressures at 
Z — 0 , with A z ~ 0. 1 , At = 0.005 , and T = 3.6 . 
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Fig. 6. Effect of number of points in x direction on harmonics at 
z = 0 , with A z = 0.1 and At = 0.005 . 


The effect of time step is shown in figure 7. A time 
step of At - 0.01 (55 steps per cycle) is sufficient to 
capture the fundamental and first harmonic but is not 
sufficient for the higher harmonics. A time step of 
At = 0.005 (110 steps per cycle) is temporally con- 
verged on this grid for all harmonics up to the third har- 
monic. 

In summary, the nondimensional grid size necessary 
for the current numerical algorithm to accurately cap- 
ture the physics of this flow is as follows: A z = 0.1 (32 
points per transverse wave) and Ajc = 0.01277 (43 
points per axial wave). The nondimensional time step 
necessary is At = 0.005 (110 steps per cycle). This 
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Fig. 7. Effect of time step on harmonics at z = 0 , with 

Ax = 0.01277 and Az = 0.1 . 


baseline grid and time step are utilized to assess the ef- 
fect of a sliding-zone interface. 

Although not shown, a two- zone grid with zone 2 held 
stationary yields a solution that is identical to the single- 
zone solution. When zone 2 slides past zone 1 at a con- 
stant velocity, the solution becomes distorted as the zone 
speed increases, as shown in figure 8. In this plot, the 
theoretical predictions and the third harmonic are not 
shown, and the computed results are displayed as lines 
rather than symbols for clarity. The zonal interface is at 
x/x = 0.5047 . Up to a speed of roughly >v Ioni , = 6 (non- 
dimensionalized by the speed of sound), the solution is 
essentially unaffected. At w zonf >12, however, the solu- 
tion deteriorates significantly. 

At a time step of At = 0.005 and a zone speed of 
w. onf = 6 , zone 2 slides past one periodic variation in 
the acoustic data (i.e., the space over which the flow 
field is periodic in the z direction) in approximately 107 
time steps. This parameter can be used as a rule-of- 
thumb indicator of the limit on the time step when a slid- 
ing zone is present. Extensive investigation by the au- 
thor (not all of which is reported here) indicates that if a 
zone is sliding past another zone and flow-field spatial 
variations in the direction of zone motion are present, 
then the time step must be set to allow approximately 80 
or more time steps for the moving zone to traverse one 
period of that variation. This rule of thumb is summa- 
rized by the following equations: 

< 18 

or 


A / 


*-pfriod 


(19) 
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Distance x/x 

Fig. 8. Effect of zone speed on harmonics at z = 0 , with 

Ax = 0.01277 , Az = 0.1 , and At = 0.005. 


where z ptriod is the spatial distance in the direction of 
zone motion over which the flow field varies by approx- 
imately one period and w d is the component of the 
acoustic-wave velocity parallel to the moving-zone in- 
terface. Because this rule of thumb is based only on 
qualitative estimates of acceptable error, it is best uti- 
lized as a guideline for establishing an initial estimate 
for the maximum time step or zone speed. Further re- 
finement to the desired tolerance can be made through 
parametric studies. 

In the case of axially-moving waves, w data = 0 . If the 
time step is fixed at A / = 0.005 , then the maximum zone 
speed according to the guideline is roughly 8. This con- 
clusion is consistent with the results in figure 8. If, on 
the other hand, the zone speed is fixed at = 30 , 
then the time step must be no larger than At* 0.0013 . 
The results obtained with At = 0.001 and w = 30 
are contrasted with those obtained with At = 0005 and 
™zone = 30 * n figures 9 and 10. The former result shows 
little distortion of the pressure contours in zone 2, 
whereas the latter shows large distortion. The oval con- 
tours should travel left to right only, with some com- 
pression, and then diminish past x/x = 1 , where the 
grid is stretched. 

This rule of thumb is independent of grid spacing in 
the direction of zone motion. Although not shown, two- 
zone results with 17, 33, and 65 grid points in the z di- 
rection yield nearly identical results when Ar and w 
are held fixed. Also, if zone 1 is clustered so that Az var- 
ies from 0.1 to 0.0001 while zone 2 retains uniform grid 
spacing, the results again follow the same rule of thumb. 
The number of x -TS multigrid subiterations currently 
employed (3) is sufficient to drive the L 2 norm of the 
equation for density down by at least 1.5 orders of mag- 
nitude for each time step. The use of six subiterations 
drives the residual down 3.5 orders for each time step, 
with no perceptible change in the solution. 





Fig. 9. Pressure contours with At = 0.001 , T = 3.6 , and 



Fig. 10. Pressure contours with At = 0.005, T = 3. 6, and 


2D Vibrating Plates 

The following 2D model problem is patterned after a 
computation by Huff. ' An infinite row of vibrating flat 
plates is simulated on the rectangular domain depicted 
in figure 11. The plates, of chord length c , are separat- 
ed by height h = h'/c' = 1 . Slip-velocity (Euler) 
boundary conditions are applied on the plates, and peri- 
odic boundary conditions are applied elsewhere on the 
upper and lower boundaries. The left and right bound- 
aries employ a characteristic boundary condition, and 
the grid is stretched in the x direction past 
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x = x/c = 5 . When used, a sliding-zone interface is 
employed at a location 0.5 chord length downstream of 
the trailing edge of the plates, which vibrate in phase 
with each other in accordance with 

z' = A'sin(o)V) (20) 

where to' is the frequency in rad/s. 


OI . ■ * -i Direction of 
Sliding interface zone 2 motion 


Zone 2 



x= 1.5 


of M - 1.2. Equation (19), with w data = ±1 (because 
the waves travel obliquely to the zone interface with a 
maximum speed equal to the speed of sound), is used to 
determine a maximum time step of At max = 0.0014 . 
Time-step studies (not shown) indicate that the time step 
necessary to adequately resolve the physics of this flow 
field is approximately At = 0.0025 , which corresponds 
to 400 time steps per period of plate oscillation. Hence, 
given the fixed zone speed of w zone = 1.2 in this case, 
the time-step constraint that is attributable to zone mo- 
tion is more restrictive than the time-step constraint 
needed to adequately resolve the flow field. 



i Plate i 


Fig. 12. Pressure contours with one zone, At - 0.001 25 , and 
7=5. 


Fig. 11. Sketch of 2D vibrating-plate case. 

If the lengths are nondimensionalized by c and time by 
c '/a'„ and the reduced frequency is defined as 


u „ 

then the nondimensional equation that defines the plate 
vibration becomes 

(M kt\ 

z = hsm\^ — - — J (22) 

Here, h is taken as 0.00004, k = 871 , and = 0.5 . 
Many acoustic waves generated by this test problem 
travel at an oblique angle in both directions relative to 
the sliding interface. At this frequency, approximately 
four complete cycles of pressure variation span the 
space between the plates. Hence, z period ~ 0.25 . A verti- 
cal uniform grid spacing of Az = 0.0078125 is used, 
which results in a spatial resolution of approximately 32 
grid points per wave. The grid spacing in the axial direc- 
tion is Ax = 0.03125 between x = -3 and x = 5 and is 
stretched past x = 5 . Although not shown, parametric 
studies indicate that this grid spacing is adequate to spa- 
tially resolve this flow field. 

Because this model problem is based loosely on the 
types of acoustic waves that may be generated by vibrat- 
ing stator vanes in an engine, computations are carried 
out with a fixed zone speed that is comparable to that of 
a typical maximum rotor tip speed. Hence, w zoflf = 1.2 
is chosen; this value corresponds to a tip Mach number 



1 Plate 1 \ 

Sliding interface 


Fig. 13. Pressure contours with two zones. At ~ 0.001 25 , 



1 Plate 1 \ 

Sliding interface 


Fig. 14. Pressure contours with two zones. At = 0.0025 , 

T = 5, and w ;onf = 1.2. 

Because of the limitation caused by zone motion, a 
time step At = 0.00125 is employed. Pressure contours 
(nondimensionalized by p ’ ja'\) for a single zone are 
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shown in figure 12. Contour levels range from 0.712 to 
0.716 in steps of 0.0002. Figures 13 and 14 show the re- 
sults of using two zones, with time steps of 
At = 0.00125 and At = 0.0025, respectively. Results 
with the smaller time step are only slightly affected by 
the sliding-zone interface, whereas results with the larg- 
er time step show greater distortion of the acoustic 
waves after they have passed through the zone interface. 
Additional details of these differences are shown in fig- 
ure 15, which is a plot of pressure levels as a function of 
x along a line at z = 0.75 . 



X 

Fig. 15. Pressures along z - 0.75 with T = 5 . 


3D. Rotor-Stator Interact ion Modes 

To assess the passage of acoustic waves through a 
sliding-zone interface in a more realistic configuration, 
theoretical rotor-stator interaction modes are derived 
based on a 16-blade 20- vane NASA Langley model,*® 
with duct radius r' = 0.1393 m, Mach number 
M = 0.6, and rotor speed 16,900 rpm. The theoretical 
levels of pressure perturbation are 

p'(r' } 0,*') = A'J m (k/r')expHa'k't'-mQ-k a 'x'] (23) 

where A' is the magnitude of the perturbation, taken as 
O.OOlpo' (Pc' is taken as 100,000 Pa), and J m is the 
Bessel function of the mth kind. One of the propagating 
modes (for a single blade passage frequency (BPF)) is 
the (-4,1) mode, where m = -4 is the circumferential 
mode number and n = 1 is the radial mode number. For 
a single BPF, k' = 87.97 m“ ! ; for the (-4,1) mode, 
= 38.17 m and k a ' = -211.37 m”*. This mode 
spins in the direction that is opposite to that of the rotor. 

To test whether this acoustic mode can be propagated 
undistorted through a patched sliding-zone interface, a 
3D time-accurate computation is performed in a cylin- 
drical duct. Because the (-4,1) mode is periodic over 
tc/2 , the grid used is a quarter of a cylinder. The grid ex- 
tends from x' = -0.5 to x' = 0 m and is divided into 
three zones. The middle zone, which extends from 
x ' = -0.2786 to x' = -0.1393 m, can be rotated. Time- 
dependent pressure perturbations are imposed at the 
downstream boundary (at jc' = 0 ) with equation (23). 


Based on grid studies (not shown), the grid that repre- 
sents the best compromise between accuracy and effi- 
ciency for this case is 273 x 33 x 33 , with Ax' = 0.00116 
m in the unstretched portion of the mesh. This grid size 
yields 33 points per wave in the circumferential direc- 
tion and approximately 25 points per wave in the axial 
direction, which is sufficient to propagate this mode 
through both downstream zones with reasonably mini- 
mal phase shift and attenuation. (Maximum pressure 
levels are computed to be approximately ±35 Pa, 
whereas the exact linear solution is ±40 Pa.) The radial 
and circumferential spacing is uniform, and the axial 
spacing is uniform in the two downstream zones (with 
121 grid points each) and is stretched in the upstream 
zone. The grid, with every fourth grid point plotted in 
the axial and circumferential directions for clarity, is 
shown in figure 16. In this figure, the middle zone is 
shown in a rotated position. 

Time-step studies (also not shown) indicate that at 
least 36-7 1 time steps per period are required to propa- 
gate this mode through both downstream zones with 
minimal attenuation and phase shift. Therefore, a nondi- 
mensional time step At = 0.001 is employed, which 
corresponds to 71 time steps per period for BPF. 


Pressure perturbation boundary 
condition applied at downstream 
boundary / 



Flow direction 


Fig. 16. 273 x 33 x 33 duct grid (points removed for clarity). 


The 3D duct acoustic modes show less sensitivity to a 
moving zone than the 2D test cases. A minimum of ap- 
proximately 40 (rather than 80) time steps is required for 
the moving zone to pass one spatially periodic variation 
in the flow field. Therefore, for rotating zones, equations 
(18) and (19) are modified to read as 

= (24) 


or 


A/ 


w period 




(25) 


A fixed moving-zone speed of b lone = 1 6, 900 rpm re- 
sults in a maximum allowable nondimensional time step 
of (with Q.. = 90° and 
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®daia = -67, 600 rpm). This time-step limitation that re- 
sults from zone motion is slightly higher than the time 
step of A / = 0.001 necessary to resolve the physics of 
the flow (regardless of zone motion). Therefore, for this 
case, the constraint due to the numerical resolution of 
the physics is more restrictive than the constraint due to 
zone motion, and the zone motion is expected to have 
little or no discernible effect on the temporally resolved 
solution. 

The acoustic waves generated by the downstream per- 
turbation propagate upstream in a spiral manner (i.e., 
like a “barber pole”). A total nondimensional time of 
T = 0.8 is sufficient to achieve a periodic solution 
through the two downstream zones. Figure 17 shows 
pressure perturbation levels as a function of x along the 
duct wall at a 0 = Constant location for a nonrotating 
zone and for a zone that is rotating at 16,900 rpm. Re- 
sults are essentially identical for both cases, which indi- 
cates, as expected, that the sliding-zone interface does 
not distort the (-^4, 1 ) mode acoustic waves with this grid 
and time step. Each computation requires approximately 
4 hr of CPU time on NASA Langley Research Center’s 
CRAY Y-MP computer. 



T = 0.8. 

Figure 1 8 shows the effect of increasing the zone ro- 
tation speed. Tripling the rotation speed results in negli- 
gible distortion, whereas quintupling the speed changes 
the character of the solution upstream of the moving- 
zone interface. With equation (24), a time step of 
At = 0.001 allows for a maximum zone speed of 53, 1 13 
rpm. This maximum speed is consistent with the results 
in the figure. Five subiterations are required for the 
higher zone speeds. 

Other propagating modes are investigated by con- 
structing a 241 x 65 x 65 grid similar to the one depicted 
in figure 16, except that Ax' = 0.000481 m and the two 
downstream zones each extend only x - 0.05 m. This 
shorter axial extent allows all modes of interest to prop- 
agate through both downstream zones in a reasonable 
computational time. This grid spacing yields between 



x\ m 

Fig. 18. Effect of zone speed on pressures along duct wall with 
At = 0.001 and T = 0.8 . 


10-65 points per wave in the circumferential direction 
and 31-82 points per wave in the axial direction, de- 
pending on the mode. The computation is performed by 
combining 10 different propagating modes at the down- 
stream boundary, with varying acoustical strengths. 
These modes are described in Table 1. Included in the 
table is the acoustical strength A' that is used in the 
computation; also included is the theoretical cutoff ratio 
of each mode. Note that these 10 modes are a represen- 
tative sample of the BPF and 2 BPF propagating modes 
in the NASA Langley Research Center model. 


Table 1 Computed Modes in 10-Mode Combined Computation 


Mode 

k ' , m 1 


V ’ m 1 

A ' , Pa 

cutoff 

(-4.D 

87.97 

38.17 

-211.37 

30 

2.88 

(-4,2) 

87.97 

66.64 

-191.80 

60 

1.65 

(-4,3) 

87.97 

91.04 

-159.55 

100 

1.21 

(-8,1) 

175.93 

69.26 

-425.84 

30 

3.18 

(-8,3) 

175.93 

127.60 

-388.83 

60 

1.72 

(-8,6) 

175.93 

200.21 

-278.67 

100 

1.10 

(12,1) 

175.93 

99.63 

-410.00 

30 

2.21 

(12,4) 

175.93 

188.41 

-306.70 

60 

1.17 

(12,5) 

175.93 

213.42 

-231.27 

100 

1.03 

(-28,1) 

175.93 

218.80 

-192.69 

150 

1.01 


Computations are performed both with and without 
middle zone rotation. For the m = -4,-8, 12 , and -28 
modes, Q period is equal to 90° , 45° , 30° , and 12.85714° , 
respectively, and Q dafa is equal to -67,600 rpm, -67,600 
rpm, +45,066 rpm, and -19,314 rpm, respectively. With 
equation (25), the time step recommended for a minimal 
effect of zone motion at Q zonf = 16,900 rpm is 
A t max « 0.0014 , 0.0007 , 0.0014 , and 0.0005 , respective- 
ly. However, in the interest of reducing the computa- 
tional time, a nondimensional time step of At = 0.001 is 
employed, which yields 71 time steps per period for sin- 
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gle BPF and 36 for 2 BPF. Because the At max from equa- 
tion (25) is exceeded for several modes, some effects 
attributable to the moving zone can be expected. Com- 
putations are run to a nondimensional time of T = 0.75 . 
The computations require 13 hr of CPU time on NASA 
Langley Research Center’s CRAY Y-MP computer or 
5.8 hr on the NAS CRAY C-90 computer. Approximate- 
ly 50 million words of memory are required. Results are 
analyzed by decomposing the computed solutions at 
various axial stations into individual modal amplitudes 
and comparing these amplitudes with the theoretical 
(undiminished) pressure amplitude. 

Results for the m = -4 , -8 , 12 , and -28 modes are 
shown in figures 19-22. Note that the numerical algo- 
rithm, with this grid and time step, has difficulty propa- 
gating those modes that have the lowest cutoff ratios, 
regardless of whether a moving zone is employed. In 
general, the lower the cutoff ratio, the more the acoustic 
wave is attenuated. Further studies need to be performed 
to determine the grid and time-step requirements for all 
modes of interest. Overall, the effects of the sliding- 
zone interface on the transmission of acoustic modes are 
very slight; the least effect occurs for the m = -4 
modes. 


240 r 

K 3 ' theory 

n 200 o no zone motion 
o Km*- 1 6.900 rpm 
Q) 

T> 160 - 
3 


Moving zone 



-0.12 -0.08 -0.04 0.00 

x ' , m 

Fig. 19. Acoustic amplitudes for m = -4 modes in 10-mode 
combined computation. 

Conclusions 

Patched sliding-zone interfaces, in combination with 
the time-accurate Euler and Navier-Stokes equations, 
are employed for aerodynamic computations of rotor- 
stator interactions (e.g., refs. 4-7). However, before 
codes with sliding interfaces can be used as a major part 
of any noise-prediction effort, the accuracy with which 
acoustic waves are passed through the sliding-zone in- 
terfaces must be assessed. Preliminary findings indicate 
that an important factor that contributes to the accuracy 
of acoustic-wave passage is the number of time steps re- 
quired for the sliding zone to move past spatial varia- 


240 r 

CO ' 

Pi 200 o 9° zone motion 

□ ^OTU ( =16 ' 900 r P m 

CD 

X> 160 - 
3 


Moving zone 



-0.12 -0.08 -0.04 0.00 

x', m 

Fig. 20. Acoustic amplitudes for m = -8 modes in 10-mode 
combined computation. 


240 p 

m ■ theory 

CL 200 o no zone motion 

□ ^on. = 16 - 900 r P m 

CD 

"D 160 - 
3 


Moving zone 



-0.12 -0.08 -0.04 0.00 

x' , m 

Fig. 21. Acoustic amplitudes for m = 12 modes in 10-mode 
combined computation. 

tions in the flow field that exist in the direction parallel 
to the sliding-zone interface (e.g., circumferential varia- 
tions for duct flows). To avoid unreasonable distortions, 
at least 40-80 time steps are required for the sliding 
zone to pass one period of the spatial variation. If the 
time step is too large, then distortions in the acoustic 
waves result. This rule of thumb has been demonstrated 
for two simple two-dimensional test problems, as well 
as for a three-dimensional-duct test problem. 
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Fig. 22. Acoustic amplitudes for m = -28 mode in 10-mode 
combined computation. 
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